A new solution to the statistics of hard elongated objects 
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We propose an analytical solution to the statistics of hard elongated objects (e.g. needles, rect- 
angles, ellipses, etc) in one dimension, by using Gibbs free energy formalism and an adequate 
approximation to the phase space. Some analytical formulation for a collective thermodynamic 
properties are found that are in agreements with the results recently obtained numerically from the 
exact model. We also verify the inverse distance law of sound pressure and notice below a certain 
pressure the expectation value of inverse distance is deviated from scaling as the inverse of distance. 
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I. INTRODUCTION 

Elasticity theory describes how a system under distor- 
tion is mechanically deformed. There are two approaches 
to quantitatively study this. The traditional approach is 
to analyze the dislocation of fluid rigid boundaries. De- 
pending on the properties of fluid (e.g. viscosity, com- 
pressibility, etc.) a variety of different cases appears, 
(see [l|]). The other approach, which is of our inter- 
est, replaces fluid with discrete objects and studies the 
short-range interaction between these objects. If the ob- 
jects are spherically-symmetric, their alignments lead to 
translational ordered/disorderedphases, @, 0, 0] • Frenkel 
et.al. in a series of papers (see [a]) developed a method 
that enables to analyze the stress and elasticity of hard 
spherically-symmetric objects. 

A natural generalization is to replace the spherically- 
symmetric objects with elongated objects. These objects 
carry a coupling between their translational and rota- 
tional degrees of freedom and display orientational or- 
dered/disordered phases; similar to liquid crystals Q. 
Recently, in a series of papers a formalism for direct 
calculation of elastic properties of hard non-spherically 
symmetric objects was proposed by Murat, Kantor, and 
Farago (see Q). They considered hard stiffness for these 
objects in order to prevent the influence of orientational 
degrees of freedom into kinetic energy. Their method was 
developed on the basis of different types of central and 
non-central short-range interactions; central potential de- 
pends only on the relative distance between particles, 
whereas noncentral potential depends on individual ob- 
ject orientations. This formalism has been so far applied 
in different problems such as the wrapping of proteins 
in DNA [f| , the ordering of complex liquids systems and 
percolation transitions [9(, and the jamming transitions, 
10]. However, this formalism is hard to be solved ana- 
lytically and a Monte-Carlo simulation should been used 
to extract its physical properties. Kantor and Kardar in 
PH proposed a self-consistency check for this formalism 
in one dimension, where instead of non-spherically sym- 



metric objects, needles are applied. The center of needles 
are on a line and the angle of each needle orientation is 
randomly chosen. They solved this model by transfer 
matrix method numerically and reported an agreement 
between the numerical and MC results. 

The purpose of our communication is to reconsider 
elongated objects in one dimension and propose an ana- 
lytical solution that, in our opinion, goes a step forward 
since it allows to obtain analytical formulation for some 
collective properties obtained so far only numerically. To 
this end, we eliminate the dependence of free energy on 
the absolute value of an angle, thus make the orientation 
completely isotropic. Interparticle distance and elasticity 
coefficients are derived into analytical formulations and 
verify the exact model results of [HI]. We generalize our 
formalism to cover different types of elongated objects 
and rep eat to derive them and verify recent results of 
(Til Il2j |. This formalism allows to evaluate other proper- 
ties of the same class, such as inverse distance between 
needles. We verify the inverse distance law of sound pres- 
sure (e.g [l4j) in high densities. As expected from apre- 
vious study on spherically-symmetric object (see [X 51 ] . in 
needles model the expectation value of inverse distance 
below a critical pressure does not scale as the inverse of 
distance. 
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FIG. 1: The phase space diagram of needles model. 

In needles model, a number of needles with centers 
located on a straight line are considered. Each needle 
has two degrees of freedom: translational x, and rota- 
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tional <fi (where —tt/2 < <fr < 7r/2). When two needles 
touch one another their orientations displays their min- 
imum distance. This distance is Ofj^+i with Xij + \ = 
sm(\(f>i ~ (f>i+i |)/ maxcos(0i, (f>i+i), and t is a length scale. 
The phase space diagram (X, </>i, 2 ) of needles model is 
depicted in Figure (p}. The two maxima at the opposite 
tips of left and right 'wings' denotes the configurations 
where two adjacent needles are bending tangentially at 
the central axis. The minima of X — denote all dif- 
ferent configurations in which two adjacent needles are 
oriented parallel, 4>i — (f>2- Every contour line in Figure 
(H|) denotes a set of degenerate configurations of a corre- 
sponding distance. The larger a distance is, the lower its 
degeneracy becomes. 

Following Takanashi's proposal (see 01) for studying 
free energy of hard objects, two adjacent needles, located 
at Xi and Xi+i interact by a hard potential U(xi, Xi + \). 
This potential is zero for larger distances than the mini- 
mal distance £X i i+ i, and infinite otherwise. 

Consider a chain of N needles compressed by the 
external pressure p. By definition the Gibbs free en- 
ergy is G = E - TS + PV. In the lack of tem- 
perature dependance, a configuration partition func- 
tion is defined as Zq = J2 e~^ G - We assume the ki- 
netic energy of particles is negligible, which means we 
study the system after all needles are at rest. Hard 
potential contributes to the partition function; Zq = 
YliLiI d 4>i J dx i exp(-f3(U(x l - x i+1 ) + p\xi - x. l+1 \)). 
Renaming Si := \xi — Xi + \\ and applying hard poten- 
tial makes up Zq = Yli=il d<t>i JpZ dsiexp(—f3psi) — 

(Pp)- n Uti I d( Pi exp(-j8j>*X(&,&+i)). Notice that 
in thermodynamic limit, open chain and closed chain of 
needles have the same free energies; (beyond this limit 
the difference are in a surface free energy that appear in 
the open chain and scales as l/N). 

In general, two systems behave likewise if their phase 
spaces are similar. This highlights acceptable approxima- 
tions that may simplify the free energy associated with 
the phase space of Figure (H|). In the following we con- 
sider a phase space in which the wrinkles of Figure ([T]) 
are flattened. This is indicated to analytically solve the 
problem and more importantly it is identical to the orig- 
inal model in the study of some collective properties. We 
can also easily generalize this formalism to also study 
some different types of elongation, other than needles. 

II. THE FORMALISM 

Consider the minimal distance between two particles i 
and i + 1 is defined by 

2 

Xii + i = - \<j>i - 4>i+i\- (I) 

The corresponding phase space diagram is depicted in 
Figure @. Comparing this with Fig. ([T]) one observes 
in the latter model the dependence of free energy on sin- 
gle angles has been vanished, therefore orientations are 



isotropic. We expect the two models behave differently 
on their ordering at high pressures; however, if we ig- 
nore the discussion of ordering and just look at proper- 
ties such as stress and elasticity coefficients the models 
behave similarly even at high compression. 




FIG. 2: Phase diagram of Eq. ©. 

Let us define the dimensionless pressure / = (3pi. 
Based on ([T]) the partition function becomes Z = 

/ 2, ,\ N 

f~ N ( / d0id0 2 e~ / - 101-021 J . In two adjacent particles 

the angles are independent and each one takes a value 
between — tt/2 and tt/2. Therefore, using Eq. (Q]), two 
needles have the distance 2x, where < x < 1. A two- 
body partition function becomes 

Zl = r'(^-)=i^. m 

Since the distance between a pair of particles is an inde- 
pendent parameter, the partition function of N particles 
is Z = Zi and the Gibbs free energy per particles be- 
comes 

(3G/N = -InZx = 2 In/- ln(l - e~ 2/ ). (3) 

The first term indicates the orientationally- 
independent free energy (point particle part) and 
the second term reflects hard interaction between the 
particles. The behavior of this free energy in different 
pressures is plotted in Figure © denoted in solid line, 
and was compared to the point particle behavior in 
dashed line. This result is in good agreement with the 
numerical results of Figure (3) in [ll| where the full 
anisotropic model was applied. The analogy of results 
highlights that the free energy of exact model of [ll| 
is unaffected if its anisotropy is eliminated. In other 
words, the major contribution into the elasticity of 
needles in thermodynamic limit is rooted at isotropic 
features. Consequently, this is not unrealistic to expect 
the parameters extracted from this free energy behave 
similar to the original model. 

Following the formalism originally proposed for 
spherically-symmetric objects in [J] and extended to the 
correlation functions of elongated objects in Q, stress 
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FIG. 3: Gibbs free energy. The solid line denotes Eq. © 
and the free energy of orientational independent particles is 
in dash line. 



1 describes a phase space almost similar to Figure ([2]) 
except that its 'wings' get bumps upwards that makes it 
look like Y-shape. In the case of a > 1 the wings are 
deformed into two valleys and the phase space diagram 
looks like a U-shape. 

In general, the partition function per particle becomes 
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Z G /N = - I e-f^ a dx 
J Jo 



— (r (I/a)- r(l/a, 2"/)), (7) 
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where T(a,x) — f r ^°t a ~ 1 e~ t dt. To get this function, we 
renamed y = f(2x) a . Using this partition function, the 
mean distance for different choices of a could be calcu- 
lated. The mean distance between the objects will be- 



and elasticity coefficient were introduced in one dimen- 
sional systems in (llj . Stress indicates mean distance 
between particles and is defined as a/ 1 — d(f3G/N)/df 
in fixed temperature T. Substituting the Gibbs energy 
from Eq. ([3]) into this definitions, interparticle distance 
becomes: 
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(4) 



The elastic coefficient is also defined by C = — da a / 9p + 
p. Re-scaling the coefficient into (3£C gives it rise to 
~~da[d] ^ f' Using the Gibbs energy in Eq. ^ the 
elastic coefficient becomes: 



C 



/ 4 cosh(2/) - (2/ 2 + / + 4) + / exp(-2/) 
2 cosh(2/) - (/a + f ) 



(5) 



The mean distance a and the elasticity coefficient C are 
plotted in Figure (g]) (see the case of a = I). The (blue) 
line of negative slope denotes a and of positive slope de- 
notes C. The dashed lines in the figure correspond to 
the mean distance and clastic coefficient of point parti- 
cles, respectively. The behavior of a and C are in good 
agreement with those extracted by numerical studies in 
Figure (4) of [HI- 



III. A GENERALIZATION INTO ELONGATED 
OBJECTS 

Lebowitz in [13j proposes a generalization into Eq. ([T]) 
of the following form: 



Xi,i + i(0i,^ 2 ) = ( -\<t>i -<k+\\ ) ( 6 ) 



Spacial cases are needles with the exponent a = I, and 
narrow ellipses with a = 2. In general, the case of a < 



a _ 1 + 1/q /i" 1 exp(-2 a /) 

f r(l/a)-r(l/a,2«/)' 



(8) 



Here we utilized the chain rule dr ( a j M ( x )) = dT(a,u) chi 

ax du ax ' 



aT(a,: 



and 

cient from a is straightforward. 



dx 



Finding the elasticity coefh- 




FIG. 4: [Color online] The distance a/l and the elastic coef- 
ficient filC versus pressure. Distances are lines with negative 
slopes and C lines are with positive slopes. The family of 
a = 0.2, q = 0.5, a = 1, a = 2, q = 5 are in order from top 
to bottom, respectively. 

C and a are plotted in the Figure for the exponents 
a =0.2, 0.5, I, 2, 5 from top to bottom, respectively. 
The mean distance lines are with negative slopes and the 
elasticity coefficient lines are with positive slops. The 
non-interacting result are depicted in dashed lines. These 
results have already been derived recently by a numerical 
analysis in 

From the agreement between the results of isotropic 
and anisotropic, one can conclude that distance and elas- 
ticity in elongated object system are highly depended 
onto the isotropic features such that even at high com- 
pressions their anisotropic effects are negligible. 



4 



IV. INVERSE DISTANCE 

Inverse distance rules the scaling of many parameters 
in physics. For instance, the inverse distance law governs 
damping of sound pressure in a medium inversely propor- 
tional to distance from sound source, [l4T |. However, in a 
granular medium (made of frictional balls) there are de- 
viation from this law, where the velocity of sound prop- 
agation scales as p v with v < 1 not p, [l5(. We shall 
analyze this deviation in the case of elongated objects. 

We calculate the expectation value of inverse distance 
directly from the free energy. The expectation value of 
inverse distance is defined: 

(a- 1 ) = - *M 

W / fdxe-f x (*) 

The derivative of left side with respect to / becomes 
| (a- 1 ) = -l + (a- 1 ) a (9) 

From Eq. ([9]) one concludes the relation (a -1 ) = 1/a is 
only verified in pressure-independent distance models. In 
other words if a depends strictly on pressure the inverse 
distance will no longer be 1/a. 

In the present model, mean distance depends on pres- 
sure. From Eq. © we should be able to directly calculate 
the mean inverse distance. Substitute Eq. (0| into Eq. 
© gives rise to the following differential equation: 

I <■-">- - 1 + (?-^t) < 10 » 

This differential equation has a solution of the form 
( -l) = / + ^l_ E e tf/\ where E 1 (z) = ^^dt 
is the exponential integral function. E\{z) is known 
to behave as a negative exponential for large / and as 
a logarithm for small /. This is standard to approxi- 
mate it by one of its bracketing function bounds (e.g. 
see Qil). In our problem, this allows to replace it with 
ie _2 ^ln(l + 1//). Consequently, the expectation value 
of inverse distance becomes: 

f Wl + 4) 

In right side of Eq. (TTTj) the first term is point parti- 
cle contribution, (and verifies the inverse distance law 
(a -1 ) ~ / in corresponding medium). The second term 
denotes the effects of elongation. In high pressures this 
correction is negligible; however, in low pressures the in- 
verse distance is gradually separated from / and gets a 



larger values. Below a certain pressure the inverse dis- 
tance approaches again to /. This transition is better 
reflected in the behavior of (a -1 ) a. This parameter in 
point particle medium is the constant 2. Elongated ob- 
jects display a deviation into it in low pressures, which is 
plotted in the inset of Figure ((5]). 
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FIG. 5: The behavior of (a ) with respect to pressure; in 
the elongated objects (solid line) and point particles (dash 
line). Inset denotes a (a -1 ) as a function of pressure, p* is 
the critical pressure. 

In conclusion, we think to have presented an accurate 
method to obtain some collective properties of elongated 
objects completely analytically. To this end, we proposed 
to eliminate the dependence of free energy onto the abso- 
lute value of an angle, and thus make the orientation com- 
pletely isotropic. This is expected to change the model 
to behave differently at high pressures, where the origi- 
nal model becomes strongly ordered; but if we disregard 
the ordering and just look at other properties (such as 
pressure or elastic constants) the differences are negli- 
gible even at high compressions. The simplicity of this 
formalism allows a better understanding of the statistics 
of elongated objects. While this formalism can receive 
applications into more realistic physical systems such as 
protein wrapping and adsorption in biology, etc., this can 
also be employed for studying elongated objects in higher 
dimensions. 

The author should thank M. Kardar and Y. Kantor for 
their valuable comments on the first draft of this note. 
He also enjoyed bright discussions with A. Nayak. This 
work was supported by NSERC Canada. 
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